Pairing Correlations on t-U-J Ladders 



> 
o 



c3 



I 

C 

o 
o 



(N 
> 

o 

m 

o 



o 
o 



X 
S3 



S. Daul 1 , D.J. Scalapino 2 and Steven R. White 3 
1 Institute for Theoretical Physics, University of California, 
Santa-Barbara CA 93106. 
2 Physics Department, University of California, 
Santa-Barbara CA 93106. 
3 Department of Physics and Astronomy, University of California, 
Irvine CA 92697. 

Pairing correlations on generalized t-U-J two-leg ladders are reported. We find that the pairing 
correlations on the usual t-U Hubbard ladder are significantly enhanced by the addition of a nearest- 
neighbor exchange interaction J. Likewise, these correlations are also enhanced for the t-J model 
when the onsite Coulomb interaction is reduced from infinity. Moreover, the pairing correlations 
are larger on a t-U-J ladder than on a t-J e g ladder in which J e ff has been adjusted so that the 
two models have the same spin gap at half-filling. This enhancement of the pairing correlations is 
associated with an increase in the pair-binding energy and the pair mobility in the t-U-J model and 
point to the importance of the charge transfer nature of the cuprate systems. 



Various ab initio quantum chemistry calculations as 
well as model Hamiltonian studies have been used to 
determine the electronic properties of Cu-oxide clusters. 
[Q-^l In particular, these calculations have provided pa- 
rameters for simpler, effective one-band Hubbard and t-J 
models which have then been used to study many-body 
correlations in larger systems. However, both the one- 
band Hubbard and the t-J models differ in an essential 
manner from the high T c cuprates which are known to 
be charge transfer insulators in their undoped state. 
Thus, the one-band Hubbard model at half-filling is char- 
acterized by a Mott-Hubbard gap which is set by U and 
in the t-J model, U is taken to infinity with the constraint 
of no double occupancy. Therefore, while Coulomb fluc- 
tuations associated with double occupancy of a site are 
controlled by U in the Hubbard model, U also determines 
the strength of the exchange coupling. In the Hubbard 
model as U increases beyond the bandwidth, J decreases 
as At 2 /U . Although J is an independent parameter in 
the t-J model, U is infinite in this model, suppressing 
charge fluctuations. 

While we believe that the basic pairing mechanism 
arises from the exchange correlations, the charge trans- 
fer nature of the cuprates can play an essential role in 
the doped systems where it allows for a more flexible ar- 
rangement between J and U than reflected in either the 
one-band Hubbard or t-J models. To explore this, we 
have carried out density-matrix renormalization group 
H (DMRG) calculations of the pairing correlations on 
two-leg t-U-J ladders. Ladders are known to provide 
model systems which exhibit various phenomena similar 
to those of the cuprates. j| In particular, when doped 
away from half-filling they are known to have power- 
law pairing correlations which have opposite, d x 2_ y 2-\ike, 
signs between the rung-rung and rung- leg correlations. 
These correlations have previously been investigated for 
both Hubbard (ujjllj and t — J models. Here we 

will study a generalized t-U-J model which includes both 



an onsite Coulomb repulsion U and a nearest neighbor 
exchange J. While both Hubbard and t-J ladders show 
pairing correlations when doped, we find that these cor- 
relations can be significantly enhanced in a model with 
both U and J. We argue that, in fact, a t-U-J model is 
appropriate for a charge-transfer material. |7]] 

The basic one-band Hubbard model is characterized by 
a one-electron nearest-neighbor hopping t and an onsite 
Coulomb interaction U . 



H=J2-t 



U 



(1) 



Here c ia creates an electron with spin a on site i and (ij) 
sums over nearest neighbor sites. As is well known, when 
U/t is large, a strong coupling expansion [Q of Eq. (|l|) 
leads to the t — J Hamiltonian 
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with J = At 2 /U and S, 5' are vectors separating nearest 
neighbor sites. Here there is an important restriction that 
no site can have two fermions. Typically in Eq.. t 
and J are treated as independent parameters and for dop- 
ing near half-filling the latter three-site term is dropped. 
Now, while these effective models both describe certain 
aspects of the cuprates system, they lack the flexibil- 
ity to describe an important feature that arises from the 
charge-transfer nature of these materials. Specifically, 
in the insulating state the one-band Hubbard model at 
large U has a Mott-Hubbard gap set by U rather than 
a charge-transfer gap set by the difference in the oxygen 
and copper sites energies. Furthermore, when U is large, 
J ~ 4t 2 /U decreases as U increases rather than saturat- 
ing at a value set by the charge-transfer gap. That is, 
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in strong coupling, the three-band Hubbard model gives 
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with i pc j the Cu (d x a_ y 2) - 0(pcr) hopping, A p( j the Cu- 
O site energy difference and U the Cu Coulomb energy. 
There are in fact further contributions to Eq. (||) from O- 
O hopping terms, as well as modifications due to O and 
Cu-0 Coulomb interactions. However, the basic point 
is that when U is large compared to the effective Cu- 
Cu hopping £p d /A p d, the exchange remains finite rather 
than going to zero. Likewise, in the t — J model, while 
J /t can be set to a physical value, one has in effect an 
infinite onsite Coulomb repulsion arising from the restric- 
tion of no double occupancy. The suppression of double 
occupancy reduces the mobility of the pairs, missing 
the physics associated with the partial occupation of the 
O sites surrounding a Cu. 

To address these limitations, we will study a t — U — J 
model in which there is a finite Coulomb interaction and 
an effective exchange term J. In the limit in which J = 0, 
this is just the one-band Hubbard model while in the 
limit U/t^> 1, this goes over to the t — J model Ea] The 
DMRG calculations reported here have been carried out 
on open ended ladders (up to 2 x 48 sites) keeping up to 
800 states, so that the maximum weight of the discarded 
density matrix eigenvalues is 10~ 6 . We first examine the 
rung-rung pair-field correlation function 



D(£) = (A l+e A\) 
for a doped (8 holes) 2 x 32 ladder. The operator 
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creates a singlet pair on the i th rung and Ai + i destroys 
it on the (i + £) th rung. A similar calculation in which a 
singlet pair is created on the i th rung and a singlet pair 
is destroyed on one of the legs at i + £ has an opposite 
sign indicating the d x 2 _ y 2-like structure of the pairing. 
Because of the finite length of the ladder, we have kept 
I < 12, with the measurements made in the central por- 
tion of the ladder, in the plots of D(£). In this region the 
effects of the open ends are negligible. 

In Fig. [j] we show the effect of adding an additional 
exchange term J to a Hubbard model with U = 6. Here 
and in the following we will measure energy in units of t. 
As seen, the addition of J clearly enhances the pairing. 
In all of the plots it is important to recognize that the 
pair has an internal structure so that Aj and Ai + g have 
only a partial overlap to the state in which a pair is added 
at the i th rung or removed from the i + £ rung, and the 
basic size of D(£) is reduced by the square of this overlap. 
As seen in Fig. 1, adding an additional exchange strongly 
enhances the pair-field correlations. 




FIG. 1. The rung-rung singlet pairing correlation function 
D(£) versus £ on a doped 2 x 32 ladder with (n) = 0.875 for 
U = 6 and various values of J. 

Similarly, in Fig. [2^,, we examined the effect of U on 
the pairing correlations of a t-U-J ladder with J = 0.25. 
For U ~S> 1, we have the usual t — J result. As U initially 
decreases, there is again a significant enhancement of the 
pairing correlations, but eventually as U decreases below 
the bandwidth, the pairing correlations are reduced. This 
is also shown in Fig. 0b, where we have plotted 



D -- 

versus U for J = 0.25. 
U ~ 6. 
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Here D reaches a maximum for 
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FIG. 2. (a) The rung-rung singlet pairing correlation func- 
tion D(£) versus f on a doped ladder with (n) — 0.875 for 
J = 0.25 and various values of U. (b) The partial singlet 
pairing correlation function sum D as a function of U for a 
doped ladder with (n) = 0.875 and J = 0.25. 

One would, of course, expect that the pairing correla- 
tions would depend on the total effective exchange inter- 
action, both the explicit "J" exchange and the additional 
exchange associated with a finite U. Thus, in the t-U- 
J model, as U initially increases, the effective exchange 
increases and then as U exceeds the bandwidth its con- 
tribution to the exchange decreases as 4t 2 /U. However, 
there is more to this than just the enhancement of the ex- 
change interaction which can be seen by comparing the 
two models. A half-filled Hubbard ladder with U = 6 
and J = 0.25 has a spin gap A s = 0.22 corresponding 
to an effective exchange Q J c g ss 2A S = 0.44. Using 
this value for the exchange in a t — J model we have cal- 
culated the pair-field correlation function D(£) in Fig. 3 
and compared it with the pair-field correlations found 
for the corresponding t — U — J model. Although both of 
these models have the same spin gap at half-filling, it is 
clear that the t — U — J ladder has significantly stronger 
pairing correlations. 



FIG. 3. Comparison of D{£) versus £ for a t-U-J model 
(solid) with [7 = 6 and J = 0.25 with a t-U model (open) 
which has the same spin gap at half-filling. 

In order to understand the reasons for this, we have 
calculated the pair-binding energy and the pair mobility 
for both these models. The pair-binding energy is 



E ph = 2E (1) - Eo(2) - Eq(0) 



(7) 



with Eo(n) the ground-state energy with n holes. We 
find -Ept, is equal to 0.34 for the t-U-J model with U = 6 
and J = 0.25. For the t-J cS ladder with J cS = 0.44, 
adjusted so that the two models have the same spin gap 
at zero doping, the pair-binding energy is 0.23. We have 
also calculated the effective hopping t c g of a hole pair 
from the dependence of 



e p (L x ) =E (2)-E (Q) 



(8) 



on the length of the ladder for ladders with L x up to 48. 
In ladders with open boundary conditions, e p (L x ) varies 
as 



e p (L x ) = e p (oo) + t cS 



(Les+IY 



(9) 



where the effective length differs from the actual ladder 
length L x because of end effects. For large enough sys- 
tems, the difference L e g — L x = SL tends to a constant 
and is considered as a fitting parameter. Fig. |J shows 
the results for the t-U-J and the t — J e g models. The 
effective hopping, given by the slope divided by 7r 2 , is 
t c g = 0.99 for the t-U-J ladder and t c g = 0.39 for the 
t — J e g ladder. 
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FIG. 4. Hole pair energy e p versus (L e g + 1) -2 for the 
t — U — J model with U = 6 and J = 0.25 (circles) and 
the corresponding t — J model with J = 0.44 (squares). The 
solid lines are least mean square fit of Eq. (^j) . 

The enhancement of the effective pair hopping which 
occurs when U is finite can be understood as arising from 
virtual states involving doubly occupied sites. An exam- 
ple of this is illustrated in Fig. 5. Here a pair of holes 
on the top rung hops to the bottom rung via a set of 
intermediate states. In this sequence, the second inter- 
mediate state, shown in the middle of the figure, has a 
doubly occupied site. In the t — U — J model this would 
not be allowed, leading to a reduction in the effective pair 
hopping. This effect not only enhances the pair-field cor- 
relations on the t—U— J ladder, but we believe also would 
act to reduce the stripe stiffness in the 2D t — J problem. 
This would favor a dx 2 -y 2 -pairing state over the striped 
state we have typically found in DMRG calculations on 
n-leg t - J ladders. §& 
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FIG. 5. Illustration of a pair transfer process involving a 
set of intermediate states. The state in the center has a dou- 
bly occupied site. This transfer process cannot occur in the 
t — J model but it does contribute to the pair hopping in the 
t — U — J model. 

Thus we conclude that the charge transfer nature of 
the cuprates can be more appropriately described using 
a t — U — J model. Furthermore, this model exhibits 
enhanced pairing correlations due to (1) an additional 
exchange coupling reflecting the exchange path in which 
there is a virtual double occupancy on the oxygen rather 
than the Cu and (2) an enhanced pair hopping allowed 
by a finite value of U which reflects the alternate paths 
for electron transfer in the charge transfer system. 
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